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Integrated power and linewidth of a propagating and a self-localized spin wave modes excited by 
spin-polarized current in an obliquely magnetized magnetic nanocontact are studied experimentally 
as functions of the angle 6 e between the external bias magnetic field and the nanocontact plane. It 
is found that the power of the propagating mode monotonically increases with # e , while the power 
of the self-localized mode has a broad maximum near 6 e — 40 deg, and exponentially vanishes 
near the critical angle e = 58 deg, at which the localized mode disappears. The linewidth of the 
propagating mode in the interval of angles 58 < e < 90 deg , where only this mode is excited, is 
adequtely described by the existing theory, while in the angular interval where both modes can exist 
the observed linewidth of both modes is substantially broadened due to the telegraph switching 
between the modes. Numetical simulations and an approximate analytical model give good semi- 
quantitative description of the observed results. 



I. INTRODUCTION 

Current-induced spin- wave dynamics^ has attracted 
great attention in recent years, both from applied and 
fundamental points of viewP"^ In terms of applica- 
tions, the possibility of realizing nanoscale microwave 
oscillators — so-called spin-torque oscillators (STOs) — is 
very appealing. From a fundamental perspective, this 
topic offers a new field of investigation in which spin- 
wave dynamics can be studied at reduced dimensions. 

In an STO, a high-density spin-polarized direct current 
(~ 10 8 A/cm 2 ) transfers part of its spin angular momen- 
tum to a ferromagnetic (FM) thin film via a spin-transfer 
torque effect P^ Under certain conditions of the applied 
magnetic field and current, a sustained precession of the 
magnetization in the ferromagnetic film can be achieved. 
Typically, a pseudo spin-valve stack is needed to both 
spin-polarize the current through one of the magnetic 
layers (referred to as the "fixed" layer), and to be able to 
read the time- varying giant magnetoresistance (GMR) 
signal arising from the magnetization precession of the 
other layer (thought of as a "free" layer). Indeed, the 
way signals are detected in these devices is as an oscil- 
lating voltage, which arises because of Ohm's law, which 
in this case can be written as V(t) = R(t)Id c - 

STOs can be realized in two geometries: nanocontacts 



(NCs), where the only nanosized feature is the nonmag- 
netic electrode through which the high density current 
is injected, and nanopillars, in which the magnetic lay- 
ers, in whole or in part, are nanosized. In this work, we 
will consider the case of NC-STOs. These devices are 
particularly interesting from a fundamental perspective, 
because they represent an ideal system in which the char- 
acteristics of current-induced spin waves can be studied. 
In fact, in such devices the spin- valve stack is patterned 
into a larger mesa with a lateral size of several microns. 
This allows for the spin waves generated at the NC to 
be damped out before reaching the edges of the mesa, 
and thus avoiding eventual reflections, which may induce 
more complex effects in the dynamics. 

A few years after the prediction of the feasibility of 
current-induced spin-wave excitations, Slonczewski pro- 
posed that the spin waves generated in a perpendicularly 
magnetized NC are propagating spin waves with wavevec- 
tor inversely proportional to the NC radius. 18 While his 
calculations described this kind of excitation well, they 
were inadequate for describing the case of in-plane mag- 
netized NCs. A theoretical solution was proposed several 
years later by two of the authors of this article,^ who 
demonstrated that experimental results could be better 
explained in terms of a self-localized spin-wave soliton — 
the so-called spin- wave bullet. 

The first studies of spin-wave dynamics in NC-STOs 



under the effect of oblique magnetic fields were first pre- 
sented in Refs. [5] and [HI However, it was with our pre- 
vious work 20 that it became possible to obtain defini- 
tive experimental confirmation of the two theoretical ap- 
proaches. We found that the localized Slavin-Tiberkevich 
mode only exists below a certain critical angle C , while 
the propagating Slonczewski mode can be excited at all 
angles 8 e . The observed excitation frequencies, thresh- 
old characteristics, and frequency agility with current 
agree with theoretical calculations and micromagnetic 
simulations ! 21 * 22 * Furthermore, the very recent observa- 
tions of spin-torque induced spin waves by means of 
micro-focused Brillouin light scattering strengthen these 
conclusions in a direct way, proving that it is the direc- 
tion of the external magnetic field that determines the 
character (localized or propagating) of the excited spin 
wav e S .™ 

Here we investigate both experimentally and by means 
of micromagnetic simulations the dependence of the 
power of the two modes on the applied magnetic field 
angle. The qualitative features of these two cases are re- 
markably different, and given the agreement between ex- 
periments and simulations, we are able to infer details 
of the magnetization dynamics. Finally, using the ana- 
lytical model that describes thermal effects in nonlinear 
auto-oscillators, we are able to describe the experimental 
angular behavior of the linewidth when only one mode 
is excited. When two modes are excited, the analytical 
model underestimates the value of the linewidth by more 
than one order of magnitude, and we introduce a new 
model in order to describe this situation. 



II. EXPERIMENTAL AND SIMULATION 
DETAILS 

In this work, the STO under investigation comprises 
a d = 40 nm circular NC on top of a pseudo spin- valve 
mesa composed of a 20 nm thick Co8oFe2o "fixed" mag- 
netic layer, a 6 nm thick Cu spacer layer, and a 4.5 nm 
thick Ni8oFe2o "free" magnetic layer. Further details of 
the material stack can be found in Ref. [25l On top of 
the STO, a coplanar waveguide (CPW) is defined to al- 
low for both injection, through a bias-T, of a current 
\Idc\ — 5-18 mA flowing from the biased free layer to 
the grounded fixed layer, and for extraction of the mi- 
crowave signal generated by the device. The external 
contact with the CPW consists of nonmagnetic ground- 
signal-ground probes (dc-40 GHz). A constant magnetic 
field noH e = 1.1 T is applied to the sample. During op- 
eration, a current \Idc\ — 5-18 mA flows from the biased 
free layer to the grounded fixed layer, injected via a bias- 
T. The microwave signal is amplified +22 dB and guided 
into a spectrum analyzer. Details of the measurement 
setup are given elsewhere.^ We calibrated our transmis- 
sion line, however, due to the wide frequency range in- 
vestigated (~ 20 GHz). We are aware that the computed 
power will not be as accurate as in the case of our re- 
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Figure 1. Power spectral density at Idc — -12.6, -14, and -16 
mA, for two different directions of the external magnetic field, 
(a) For e = 30 deg, two distinct peaks can be identified: a 
red-shifted lower frequency peak, identified as a localized soli- 
ton mode, and a blue-shifted higher frequency peak, identified 
as the propagating mode, (b) At e = 70 deg, only the prop- 
agating mode is observed. 



cent work, 26 where the frequency range was much smaller 
(~ 500 MHz). We therefore expect a 3 dB error in power 
(i.e. a factor 2) to be unavoidable in the present mea- 
surements. Nevertheless, given the variation of the STO 
power over several order of magnitudes as the magnetic 
field angle is varied, the measurements are meaningful 
and sufficiently accurate in this context. 

In order to compare our experimental results with mi- 
cromagnetic simulations, we used the data from Ref. [22j 
The sample geometry in that case was very similar, al- 
though the applied field was slightly different. Since we 
are only interested in a qualitative understanding of the 
results obtained by using this experimental setup, a new 
full-scale ad hoc numerical study was not performed. De- 
tails of the micromagnetic framework can be found in 
Ref. HI 



III. RESULTS 

Fig. [I] summarizes the results of our previous work. 20 
In particular, it shows that for a certain angle 6 e < 6 C 
(0 C « 55 deg), two spin- wave modes with opposite fre- 
quency agility can coexist. Only one mode exists for 
e > C . Comparison with theory and micromagnetic 
simulations allowed us to identify the mode at lower fre- 
quency as a localized mode of solitonic character (a spin- 
wave "bullet"), and the high frequency mode as a prop- 
agating modeP^ 

Fig. [2|a) shows the measured integrated power of the 
two modes as a function of applied magnetic field an- 
gle. The integrated power is calculated by averaging the 
power at each current above the threshold current.^ The 
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Figure 2. (Color online) (a) Averaged experimentally mea- 
sured integrated power of both modes as a function of angle 
of applied magnetic field. The asymmetric error bars at each 
data point represent the upper and lower range of the data set, 
which is made up of the integrated power at different current 
values, (b) Non-linear power of both modes computed from 
micromagnetic simulations using Eq. |l]). 



threshold current is obtained by the procedure illustrated 
in our previous work™ Error bars indicate the range 
over which the power varies at different bias currents. 
The propagating mode shows an approximately mono- 
tonic increase in emitted power as out-of-plane angles 
are approached. The localized mode, on the other hand, 
shows a monotonic increase in power up to an angle of 
e « 30 deg, followed by a plateau up to e ~ 45 deg, and 
then a rapid decrease down to the critical angle, above 
which the mode is no longer excited. 

For micromagnetic simulations, the power was com- 
puted employing the definition of nonlinear oscillation 
power developed in Ref. (27[ and by then performing an 
analogous averaging over the current values. In Fig.[2|b), 
we show the power computed by means of this technique, 
which will be discussed in more detail in the next section. 
It should be noticed that the power of the propagating 
mode increases monotonically with the applied field an- 
gle, while the power of the localized mode increases up 
to the maximum at e « 50 deg, and then undergoes a 
decrease down to the critical angle. 

Fig. [3] shows the measured linewidth as a function of 
magnetic field angle e . The asymmetric error bars on 
each data point represent the upper and lower standard 
deviations of the data set, which is made up of the in- 
tegrated power at the various current values. The local- 
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Figure 3. (Color online) Linewidth of both modes as a function 
of angle of applied magnetic field. The continuous line is the 
linewidth calculated using the analytical model of Ref. I281 



ized mode has its minimum linewidth at e ~ 30 deg, 
while for the propagating mode the linewidth minimum 
is at e ~ 60 deg. When approaching this angle at which 
the localized mode disappears, the linewidth drops very 
sharply, by almost two orders of magnitude, as the angle 
is varied by « 2 deg. 



IV. DISCUSSION 

A. Integrated power 

In order to describe the angular dependence of the 
integrated power, we use a mathematical model based 
on the evaluation of the nonlinear oscillation power™ 
Such a quantity, which represents a proper generaliza- 
tion of oscillation power to the case of spatially nonuni- 
form and noncircular magnetization precession, is indeed 
found to be proporti onal to the integrated power com- 
puted experiment allyP^II This model requires knowl- 
edge of the closed trajectory of the magnetization vector 
during precessional motion. By normalizing this vector 
w.r.t. the saturation magnetization, the dynamics occurs 
on the unit sphere. The corresponding curve splits the 
surface of the unit sphere into two parts: a larger part 
Sl, and a smaller part Ss- The nonlinear power Pnl is 
defined as the sum over all computational cells N of the 
areas of the smaller surfaces Ss : 



NL 



\rtff dS 



4irN 



(1) 



It should be recalled that, since the simulations per- 
formed in Ref. [22] were based on the realistic assumption 
that the spin-transfer torque only dynamically affects the 
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Figure 4. (Color online) Precession orbit of the magnetization 
for the localized mode (top row) and the propagating mode 
(bottom row), at different applied field directions, indicated 
by the solid arrow. The magnitude of the field is the same in 
all cases. 



magnetization configuration of the thinner "free layer" , 
the formulation in Eq. (HI) will therefore not account for 
the trajectories described by the magnetization vector of 
the thicker "fixed layer". 

The results presented in Fig. [2] show that the qual- 
itatively different experimental behaviors [Fig. [2fa)] of 
the power of the two modes as a function of the applied 
field angle e (a monotonic increase for the propagating 
mode, a relative maximum for the localized mode) can 
be well reproduced by micromagnetic simulations [Fig. 
[2|b)]. Quantitative agreements are difficult to obtain for 
several reasons. Firstly, the discrepancies between the 
nominal and realistic parameters values are hard to esti- 
mate. Secondly, some relevant contributions arising from 
finite temperature effects, Oersted fields, and nonlinear 
damping were not considered in the simulations. Recent 
numerical investigations^ indeed show that a satisfying 
quantitative agreement between micromagnetic simula- 
tions and laboratory experiments can be only achieved in 
the presence of a large Oersted field. In this framework, 
it was demonstrated that the inhomogeneities created by 
this field contribution drastically change the properties 
exhibited by both the propagating modes (the frequency 
tunability, for example, can even change sign 20 -*) and lo- 
calized modes (the threshold current increases 29 ). One 
may therefore ask whether simulations that do not in- 
clude the Oersted field are reliable at all, even if only for 
qualitative comparison. 

For this reason, we performed simulations including 
the Oersted field at two selected angles 9 e = 25 deg and 
e = 70 deg. At e = 25 deg, we again found the same 
qualitative behavior as in our recent work, 20 namely that 
the two spin-wave modes are excited alternately. When 
we compared the precession orbits with and without the 
inclusion of the Oersted field, we observed that the effect 



of the Oersted field is one of a slight perturbation. At 
6 e = 70 deg, only one mode can be excited and, again, the 
precession orbit is only slightly changed, which in turn 
leads to a slight change in the computed power. This is 
shown in Fig. [2] by open square and filled star symbols 
for the spin- wave bullet and the propagating mode, re- 
spectively. For the spin- wave bullet there is no noticeable 
change in the computed power, while for the propagat- 
ing mode the computed power is slightly lower than the 
case without the Oersted field. The conclusion is there- 
fore that the Oersted field, while having a fundamental 
role in triggering the switching between the two spin- 
wave modes, does not perturb in a critical way the tra- 
jectory of the magnetization when a particular spin- wave 
mode is excited. The reason for the discrepancy between 
the experimental and the computed power could there- 
fore lie in (besides the aforementioned realistic parameter 
estimation and finite temperature effects) the different 
dwell times of the two spin-wave modes at different an- 
gles. However, it is beyond the scope of this article to 
investigate such effects. 

From the results of the micromagnetic simulations, we 
are also able to infer the details of the magnetization 
dynamics. For the propagating mode, the magnetization 
vector describes a surface on the unit sphere whose area 
increases monotonically with the increase of the out-of- 
plane bias field angle e , as shown in Fig. [21(a) -(c) . This 
leads to the increase in power with 6 e observed in Fig. 
|2J Such a result can be understood considering that the 
demagnetizing field in the thin film is larger at larger 
e , acting so as to open up the precession cone and to 
enlarge the corresponding precession surface. This is also 
confirmed by the fact that, at a given bias current, the 
precession frequency decreases when e increases! 20 * 30 * 

Analogous conclusions could be drawn by using a 
previously-proposed simple model for the integrated 
power, based on the assumption that the precession orbit 
is circular. 5 Under this condition, the power P can be re- 
lated to the maximum GMR signal AR max , to the angle 
between the precessional axis direction in the "free" layer 
and the magnetization in the "fixed" layer 7p, and to the 
precession angle /?, through the following relationship: 



P ~ AR max sin 7p sin f3. 



(2) 



By performing a similar numerical investigation as a 
function of ei one can also conclude in this case that the 
main effect of the increase in the out-of-plane angle is to 
increase the precessional cone by means of an increase in 
the angle /?, which in turn determines the increase in the 
power P. 

However, it should be pointed out that the linear ap- 
proach employed in Eq. (J2| is not suitable for capturing 
those details of the magnetization dynamics which are 
associated with the excitation of nonlinear localized bul- 
let modes. In fact, the closed trajectories described on 
the unit sphere are, in this case, not circular in the whole 
range of e . Elliptical clamshell-like curves are described 
instead. In addition, for large supercriticality, the applied 



current does not yield a further increase of the preces- 
sion cone, but rather induces distortions in the shape of 
the surface. Since the computed output power necessar- 
ily takes into account all these mechanisms, the nonlinear 
oscillation power tool of Eq. ([!]) appears to be more ap- 
propriate. 

In this case, the trend observed experimentally for the 
integrated power can be related to the expansion and 
contraction of the precession surface described at dif- 
ferent out-of-plane angles. In particular, for field direc- 
tion e w deg, the magnetization precession takes place 
about an effective field which is antiparallel to that di- 
rection, and the corresponding surface turns out to be 
relatively small: see Fig. ^[c). The magnetization trajec- 
tory undergoes an enlargement as the out-of-plane angle 
increases up to e « 50 deg. In such a configuration, the 
unit vector describes an elliptical quasi-diametric path 
along the sphere (Sl ~ Ss), which yields a maximum in 
the nonlinear oscillation power: see Fig. |2|d). Any fur- 
ther increase of e induces a reduction in the precession 
surface w.r.t. the previous quasi-diametric path. Indeed, 
the nonlinear power exhibits a rapid decrease down to 
the critical angle e « 60 deg: see Fig. [2je). At this field 
orientation, the magnetization trajectory associated with 
localized modes would become null, making excitation of 
these modes inconceivable. 

Above the critical angle, all the power injected by the 
current to the system is provided to the propagating 
mode. 



B. Linewidth and peak power 

In order to properly understand linewidth broadening 
in STOs one has to include thermal effects, which were 
neglected in our micromagnetic simulations (where T = 
K). Those simulations will therefore not be used in the 
following. Instead, we rely on two different models which 
apply to cases where either one or two spin-wave modes 
can be excited. In the first case, we will use the nonlinear 
auto-oscillator theory developed in Ref. [28j while for the 
second case we will introduce a new model in which the 
concept of partial coherence of the oscillator is used. 



Q > is a phenomenological coefficient which quanti- 
fies the nonlinearity of the damping! 31 * 32 ! The explicit ex- 
pressions for N can be found in Ref. [32] j p is the angle 
between the stationary magnetization of the "free" and 
"fixed" layers, calculated by solving the magnetostatic 
boundary conditions, e is the modulus of the electron 
charge, and h is the reduced Planck constant. 

Compared to a conventional auto-oscillator, the extra 
term involving the N/T e s ratio (whose value can vary 
from -5 to 30 with our material parameters) in Eq. (|3| is 
the one that governs linewidth broadening in nonlinear 
auto-oscillators. The angular dependence of the linewidth 
is also predominantly in the same N/T e s ratio, since To 
and j3 are always positive and typically vary by 50% over 
the entire angular range. One therefore expects a mini- 
mum of the linewidth at the angle where the N/T e ^ ratio 
is minimum. N is a monotonic function of the applied 
field angle in the range 0° < e < 90°, and it has a nega- 
tive value at e =0°, and is positive at e = 90°. There- 
fore, there exists an angle at which N = and where A/ 
is minimum. This angle coincides with the critical angle 
C , since the localized mode can be excited only when 
7V<0. 

In Fig.[3j we plot the theoretically calculated linewidth 
as a continuous line. The numerical constants were set 
as 7 = 1.76 x 10 11 Hz/T, fi M = 0.8 T, a G = 0.01, 
V = ttRI • t, with R c = 20 nm and t = 4.5 nm, 
which are the nominal values for our sample. We also set 
Q = (2ujm/^o) — 1 according to the theory, and £ = 1.5, 
which is rather an intermediate value for the currents 
considered here. The two free parameters are the effec- 
tive volume V^ff, which is assumed to be 10 times the vol- 
ume of a cylinder with the diameter of the nanocontact 
and height equal to the thickness of the "free" magnetic 
layer, while e = 0.2 is chosen as reasonable approxima- 
tion to the spin-polarization coefficient. From Fig. [3J it 
is evident that the analytical theory correctly predicts 
the occurrence of the linewidth minimum at the criti- 
cal angle C . The theory can also describe the data for 
e > 6 ' c , where only one spin- wave mode is excited and 
the linewidth increases at larger angles. However, the an- 
alytical model underestimates the value of the linewidth 
by almost two orders of magnitude at angles e < C , 
where both modes are excited. 



Nonlinear auto- oscillator theory 



2. Partial coherence model 



According to the model in Ref. 28, the linewidth of a 
nonlinear auto-oscillator can be expressed as 



Af= To(kBT\ 
J 2ir\/3P J 



N 

r e ff 



(3) 



where T = T G (1 + QP ), = fi M uJoV^/j, r eff = 
(a G uj Q + a I), and a = (ejh/2eM V e ^)cosj p . T G 
is the Gilbert damping rate (ferromagnetic resonance 
linewidth), Po = (( - 1)/(C + Q), where ( = I/I th , and 



In order to explain the origin of the generation 
linewidth (500 MHz) in the bimodal regime of generation, 
we start from two intuitive considerations, and from from 
the results of the micromagnetic simulations performed 
in Ref. H3 

The first consideration is that we observed hopping 
between the two modes. Assuming that this process is 
random (a usual assumption for hopping), the linewidth 
of each mode would be equal to the inverse of the mode 
lifetime r w 0.3 ns, A/ - 1/r w 3.5 GHz. Clearly, this 
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Figure 5. Time profile of individual mode oscillations in the bi- 
modal regime of STNO generation (schematically). The time 
profile consists of a series of "bursts" of approximately the 
same shape at the moments t n . The oscillation phase at each 
burst is 6 n . 



interpretation overestimates the linewidth by almost an 
order of magnitude. 

The second consideration is that at K, the hopping 
between the two modes has a well-defined pattern. That 
is, the hopping is not random, but is a deterministic pe- 
riodic process and therefore does not lead to linewidth 
broadening. The spectrum at K should consist of a 
series of (5-peaks at multiplies of the hopping frequency 
1/T. At finite temperatures, each spectral peak will be 
broadened by thermal fluctuations by about the same 
amount as the linewidth in a single-mode regime. There- 
fore, the generation linewidth can be estimated using the 
single- mode theory, Af « 20 MHz. Clearly, this approach 
underestimates the linewidth by more than one order of 
magnitude. 

The failure of both the "completely random" and 
"completely deterministic" pictures of mode hopping 
suggests that, in reality the mode hopping is only par- 
tially coherent That is, the process of persistent hopping 
happens even at K and is deterministic. However, the 
power of each mode between "bursts" reduces by sev- 
eral orders of magnitude (see Fig. I5J. When the energy 
of the mode is very low, it becomes much more sensitive 
to thermal fluctuations (the influence of thermal fluctu- 
ations is proportional to ksT/E, where E is the mode 
energy). As a result, the phase <j) n of the next burst will 
not be exactly the same as the phase of the previous 
burst n _i, but will have a significant random compo- 
nent, (j) n = (j) n -\ + 5(j) ni collected during the low-power 
stage of the evolution of the mode. In principle, the po- 
sitions t n of the bursts' maxima will also be randomized, 
but we will neglect this effect for simplicity — it has the 
same qualitative influence on the linewidth as does the 
randomness of the phases <j) n (and there are reasons to 
believe that the randomness in t n will be smaller than 
the randomness in cb n ). Thus, here we will assume exact 
amplitude periodicity: t n = nT. 

Under these assumptions, the time profile of the slow 



(without trivial dependence e~ ZUJot , where uq is the av- 
erage frequency of oscillations) complex mode amplitude 
can be written as 



c(t) = ^2a(t-nT)e i4 



(4) 



where a(t) describes the amplitude profile of one burst, 
and ct> n is the random phase of the n-th burst. We will 
assume that the phase difference (cj> n — (j) n -i) of two con- 
secutive bursts is distributed with probability 
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Here Acj) is the average phase variance between two bursts 
(phenomenological parameter). The probability distribu- 
tion Eq. ([5| corresponds to a "random walk" of the phase, 
and is valid for Acj) < 2tt. 

Using Eq. Q, one can find the Fourier image of the 
complex amplitude, 



c -= J2 



a u e 



inujT-\-i(f) n 



(6) 



where a w is the Fourier image of the amplitude a(t) and 
the energy spectrum of the oscillations is (note that since 
we are using a slow amplitude, the spectrum is measured 
from the central frequency ujq) 



S(lj) = S (oj) 



E 



R k e 



ikujT 



(7) 



Here S (uj) = \a u \ 2 /T, R k = (exp[z(</> n+fc - <j> n )\), with 
the angular brackets denoting averaging over the statis- 
tics of phase fluctuations. For the random walk process 
Eq. ([5|, the phase correlator can be easily calculated: 



ft„=Cxp(-M A(/) 2 



(8) 



Then, the summation in Eq. ¥A can be performed ana- 
lytically, yielding the final result for the oscillation spec- 
trum: 



S(uj) = coth(A^ 2 / 4 ) Sq(u) ct^uj) 



where 



o-<t>(uj) = 



cosh(A^ 2 /2) 



cosh(A</> 2 /2) - cos(cjT) 



(9) 



(10) 



The oscillation spectrum S(u) is a product of two 
terms — the spectrum of the individual burst Sq(w), and 
the effective spectrum of phase fluctuations 0-^(00). The 
spectral width of the burst spectrum Sq(uj) is of the order 
of 1/r « 3.5 GHz. The spectrum a^ui) is shown in Fig. [6] 
for several values of the phase variance Acj). One can see 
that for small Acj) the width of cj^ is much smaller than 
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Figure 6. Effective spectrum of phase fluctuations cr^uS) for 
several values of the phase variance A<j>. 



1/T, whereas for large Acf) it is larger than 1/T. Two 
regimes of generation can be distinguished. 

In the first (decoherent) regime, A<j) ^> 1, different 
bursts are not correlated, and the spectrum of oscilla- 
tions is determined by the width of the single-burst spec- 
trum So(uj). This is the standard hopping regime, which 
is clearly not realized in our case. 

In the second (partially coherent) regime, Ac/) < 1, 
the phases of consecutive bursts are partially correlated, 
and the generation linewidth is determined by the width 
of the phase spectrum a^(uj). Using a series expansion 
of cr^u), one can derive an approximate expression for 
FWHM linewidth, valid in the limit where A</> <C 1: 



A/: 



2ttT 



(11) 



Using this expression and the experimental value Af = 
500 MHz, one can estimate A<fi = 1.33 = 76° (the nu- 
merical solution, without series expansion, gives the close 
value Acj) = 1.29 = 74°). 



V. CONCLUSIONS 

We measured the power emitted by an NC spin-torque 
oscillator as a function of out-of-plane applied magnetic 
field direction. As the angle varies, two distinct spin- 
wave modes are excited: the Slonczewski propagating 



mode and the Slavin-Tiberkevich self-localized mode (the 
spin- wave bullet). The power of the propagating mode 
increases monotonically as out-of-plane angles are ap- 
proached, while the power of the localized mode reaches 
a maximum, and then decreases towards a critical angle, 
above which the localized mode no longer exists. Both 
facts are qualitatively reproduced by micromagnetic sim- 
ulations. 

From the micromagnetically simulated nonlinear oscil- 
lation power, we were able to infer details of the pre- 
cession trajectory described by the normalized magne- 
tization vector on the unit sphere. In particular, we 
found that the trajectories followed by localized modes 
increase, on average, up to an out-of-plane angle e ~ 50 
deg (where an almost diametric path is described), and 
then undergo a rapid decrease down to the critical angle 
C « 60 deg (where the magnetization trajectory shrinks 
towards a single point). Above this critical angle, the ex- 
citation of bullet modes is therefore prevented, and all 
the energy provided by the current to the system only 
supports the excitation of propagating modes. 

Using the nonlinear auto-oscillator model that includes 
the effect of temperature in the calculation of the gener- 
ation linewidth, we were able to quantitatively describe 
the experimental observations when only one spin-wave 
mode is excited. When two modes are excited, the model 
underestimates the experimental linewidth by approxi- 
mately one order of magnitude. Intermediate value of the 
generation linewidth Af = 500 MHz, when both spin- 
wave modes are excited, could be explained by assum- 
ing partial phase coherence between consecutive bursts. 
The estimated value of the phase variance, A<p c± 75°, 
seems reasonable taking into account the substantial re- 
duction in the power of the mode between the bursts. 
Since the magnetization dynamics in the bimodal regime 
is extremely complex, it is unlikely that the value of A(j) 
can be calculated analytically. 
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